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FOREWORD 


The work described in this report was performed by the United Technologies 
Research Center for the NASA-Lewis Research Center under Contract NAS3-23273. The 
program manager was Dr. Anthony J. Dennis and the principal investigator was Dr. 

Brice N. Cassenti. The NASA-Lewis technical project manager was Dr. Robert L. Thompson. 
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1.0 PROGRAM OUTLINE 
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This is the Final Report for the Research and Development Program for the 
development of Advanced Time- Temperature Dependent Constitutive Relationships. 

This program is being conducted by United Technologies Research Center for the 
NASA-Lewis Research Center under Contract No. NAS3-23273. Volume I contains the 
theoretical discussion and Volume II is the Programming Manual. 

The objective is to improve the accuracy and efficiency of the nonlinear 
constitutive theories (computer models) , developed and evaluated in Contract 
NAS3-22055. The program is divided into six tasks. The first task concerns 
development of a theoretical formulation to allow temperature dependence of the 
equilibrium stress state variable during elastic excursions in a thermomechanical 
hysteresis loop. In the second task, a new term is incorporated into three candi- 
date theories to accommodate an inelastic instantaneous response. In the third task, 
the combined effects of Tasks I and II are investigated. In the fourth task, 
various techniques are incorporated into the subroutine HYPELA of the MARC*program 
to speed up computation. The fifth task involves the verification, documentation 
and delivery of the software program. The sixth task is reporting. 

The Task I objective is to develop a theoretical formulation which allows for a 
change in the equilibrium stress state variable with temperature during elastic 
excursions in thermomechanical hysteresis loops. This formulation is incorporated 
into the three nonlinear constitutive theories evaluated in the previous contract 
NAS3-22055 and is encoded into the MARC* three-dimensional nonlinear finite 
element program. The comparison of these theories is completed by analyzing the 
response using three thermomechanical hysteresis loop tests performed on Hastelloy-X 
material. 

The Task II objective is to incorporate into the candidate theory selected in 
Task I, a new term to accommodate an inelastic instantaneous response. For compara- 
tive purposes, the three thermomechanical "faithful cycle" hysteresis loops analyzed 
in Task I are repeated and analyzed in this task. 

The Task III objective is to evaluate the theories developed to determine the 
combined improvement in predicting stresses and strains. The analytical results of 
this task are to be compared with those of Tasks I and II and Contract NAS3-22055 
using the same thermomechanical hysteresis loops and geometry models used in Task I. 

The Task IV objective is to improve the computational efficiency by: 


* 

MARC is a three-dimensional nonlinear finite element program available from 
MARC Analysis Research Corporation, Palo Alto, California. 
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(1) incorporating an array containing the state variables in the viscoplastic 
theories into subroutine HYPELA of the MARC program to speed up computa- 
tion time, 

(2) further reducing computation time by modifying the MARC program so that 
the stiffness matrix is assembled once only during a calculation, with no 
further reassembly during succeeding increments. Furthermore, more 
efficient programming of the computer models developed in Tasks 1 and II 
and Contract NAS 3-22055 is completed, and 

(3) performing test cases of sufficient complexity to adequately demonstrate 
the improvements in computational efficiency. 

Task V concerns verification, documentation and delivery of the software program. 

Task VI consists of bimonthly technical and monthly financial reports together 
with a final technical report. 
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2.0 SUMMARY OF RESULTS 


During the course of this Investigation, several significant conclusions 
were drawn. They are: 

1. The variation of material parameters with temperature must be carefully 
considered. This variation is best illustrated by examining the variation of the 
elastic moduli with temperature. If instantaneous, or tangent, moduli are used, a 
perfect elastic reaction where the material undergoes changes in temperature cannot 
be simulated. If secant moduli are used, then elastic nonisothermal reactions can 
be simulated. Similar considerations must be examined for other material parameters 
in any constitutive model. The use of secant material parameters in time dependent 
reactions introduces terms containing rate of change of temperature. For eximple, 
during nonisothermal elastic reactions the rate of change of stress will include 
strain times , the derivative of modulus with respect to temperature times the 
derivative of temperature with respect to time. 

2. Including rate of change of temperature terms in the equilibrium stress 
formulation can have a significant effect on the predicted stress ratcheting during 
cyclic thermomechanical fatigue loading. The thermal ratcheting was also found to 
be sensitive to the material parameters used. 

3. An instantaneous (i.e., plastic) response term has been introduced which 
can be shown to reduce to classical perfect plasticity. This term allows more 
accurate reproduction of low temperature material response. 

4. The importance of accurately evaluating the material parameters is illus- 
trated by comparing the poor comparisons at the end of Task I with the improved 
results at the end of Task III where more accurate material parameters were intro- 
duced . 


5. A sixty percent reduction in CPU time can be attained by integrating the 
constitutive equations using a variable time step when compared to a fixed time 
step. 


Conclusions 1 and 2 are discussed in the Task I description. Conclusion 3 
is discussed in the Task II description. Conclusion 4 is illustrated by comparing 
the Task I and Task III results (in Task III, the material parameters were changed), 
and Task IV illustrates Conclusion 5. 
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3.0 TASK I - TECHNICAL DISCUSSION: 
RATE OF CHANGE OF TEMPERATURE EFFECTS 


The Task I objective is to develop a theoretical formulation which allows for a 
change in the equilibrium stress state variable with temperature during elastic 
excursions in thermomechanical hysteresis loops. This formulation is incorporated 
into the three nonlinear constitutive theories evaluated in the previous contract 
NAS3-22055 and encoded into the MARC three-dimensional nonlinear finite element 
program. The three theories considered are discussed in detail in Ref. 1, along 
with several others, and include Walker's Theory, Miller's Theory and Kreig, 
Swearengen and Rhode's Theory. The predicted material response from each of the 
theories, using three thermomechanical hysteresis loop tests, is compared with 
results obtained on Hastelloy-X material. 

This section consists of several subsections. Section 3.1 presents the three 
thermomechanical hysteresis loops modeled under this contract. The next section 
(3.2) discusses the variation in elastic moduli with temperature while the third 
section (3.3) introduces the three constitutive theories considered and the modifi- 
cations that were made to include rate of change of temperature terms. The next 
three sections (3.4, 3.5 and 3.6) compare the predicted results for each of the three 
theories, with and without the rate of change of temperature terms. 


3.1 Thermomechanical Fatigue Loops 

Three thermomechanical fatigue (TMF) loops are used in the evaluation of the 
modified viscoplastic constitutive formulations to predict the response of Hastelloy-X 
combustor liners. These loops, with their imposed loadings, and their experimentally 
measured hysteresis loops are presented in Figs. 1-3. In the first and second loops, 
Figs. 1 and 2, the temperature varies sinusoidally with time, from 427 to 871 C, 
with a period of one minute. In the third loop, Fig. 3, the temperature again 
varies sinusoidally for one minute, going from 954 C to 504 C and back to 954 C. 

The temperature is then held for 40 seconds at 954 C which completes one cycle; 
the strain varies according to Fig. 3. 

The thermomechanical loops used in this investigation represent a gradual 
increase in complexity. Figure 1 is a closed loop with equal maximum positive and 
negative strains. Figure 2 is an open loop but retains the equal maximum positive 
and negative strains while Fig. 3 is open with unequal maximum and minimum strains. 
These three cycles taken together will allow more accurate determinations of response 
of gas turbine engine hot section components. 
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3.2 Rate of Change of Temperature Effects 
in Thermoelasticity 

The three TMF loops illustrated in Figs. 1-3 include not only changes in 
mechanical loading (i.e., imposed strain) but also large changes in temperature. 
To illustrate the effects of rate of change of temperature, consider an elastic 
response only. Three forms of the constitutive equations will be considered . •, 
dimensional form. These three forms are compared to show that if perfect r ’ ?_s- 
ticity is to be represented then race of change of temperature terms will appear. 
The three forms include: 

(1) a tangent moduli formation 


dr = E(«-aT' 

(2) a secant moduli formation 

<r=-E[«-a(T-T 0 )] 

and 

(3) a mixed formulation 

where E * E(c,T) is Young’s Modulus 

a ■ a(e,T) is the coefficient of thermal expansion 
(’)- is the derivative with respect to time 

o is the stress 

c is the strain 

T is the temperature, and 

T Q is the stress free temperature. 

Equation (1) is the form presented for Miller's Theory and Krieg, Swearengen and 
Rhode's Theory in Ref. 1 for the elastic response, while Eq. (3) is the form pre- 
sented in Ref. 1 for Walker’s Theory. The MARC code assumes a form as in Eq. (3) 


( 1 ) 


( 2 ) 
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For an elastic response the material, by definition responds in a reversible 
manner, that is the stress is a unique function of strain and temperature, or 


o- = 0-(«,T) 


( 4 ) 


Differentiating Eq. (4) with respect to time and comparing to Eq. (1) yields 


is 

d* 


= E 


and 


d<r 

dl 


= -a E 


( 5 ) 


( 6 ) 


The mixed partials, Eqs. (5) and (6) must be equal. Thus, 


dE 

dT 


+ 


d(aE) ^ 

- s — *0 


( 7 ) 


If the material is linear in strain in the elastic region, E is not a function of 
strain and Eq. (7) becomes 


j_ dE + da 
E dT dt 


*0 


( 8 ) 


This is a severe requirement for the coefficient of thermal expansion to satisfy 
and is based on a uniaxial stress state only. The more general case of a multi- 
axial stress state presents even more severe requirements. This case is examined 
in Appendix 1. 


For the multiaxial stress states, Eq. (7) becomes, from Appendix 1, 


d(«K) 8jj dr 
a«ij 3 dT 


*0 


(9) 
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where K 


E 

3(l-2v) 


Is the bulk modulus, and 


6 


ij 


1 i - j 

0 1 + j 


Equation (9) is equivalent to Eq. (7), but an even more severe requirement 
on the shear modulus, G, is shown to hold in Appendix 1. 


d& 

dJ 


0 


( 10 ) 


Equation (10), requires that the shear modulus is not a function of temperature if the 
stress-strain law is elastic (i.e., if the stress-strain law is path independent). 
Therefore, it is not possible to represent an elastic response using temperature 
dependent elastic moduli defined as in Eq. (1). 

The secant moduli formulation, Eq. (2), obviously satisfies the requirement 
represented in Eq. (4). Assuming the modulus and the coefficient of thermal expan- 
sion do not depend on strain 


The rate of change of stress in Eq. 

cr = EU-aT) 


E = E(T) and 
a = a( T) 

(2) becomes 



d(aE) l^ 
dT ■» 


( 11 ) 


Tne first term is the same as that in Eq. (1). The second term is proportional 
to the ’-ate of change of temperature and illustrates that rate of change of tem- 
perature terms can enter into a constitutive formulation in addition to the effects 
of thermal expansion. 

The mixed formulation, Eq. (3), also places constraints on the material formu- 
lation. Taking the derivative of Eq. (3) with respect to time 




( 1 ?) 
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From Eq. (3) , 


and substituting into Eq. (12) 





iiiit 
E dT J T 


For an elastic response, the stress is a single valued function of the 
and temperature, therefore. 


a - cr («,T ) 


and 


da . , da . 

r 37 ,+ 5t t 


Comparing Eqs. (14) and (16) 


da _ _ , a_ da 
d * -t E d< 


and 


d<r ad* 

ar -«*Tar 


The mixed second partial derivatives with respect to strain and temperature 
be equal, which, from Eqs. (17) and (18) becomes. 


dE I a^aE _ a(aE) t I do- dE 
dT E dT dt dt E dc dT 


(13) 

(14) 
strain 

(15) 

(16) 

(17' 

(18) 

must 

(19) 
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Substituting Eqs. (17) and (18) into Eq. (19) results in 


da 

de 


0 


( 20 ) 


Appendix 2 shows that for the case of multiaxial loading, Eq. (20) becomes 


■£*-0 


( 21 ) 


and there are no other constraints on the material constants. 
Equations (]) and (2) can be recast in integral form as 


and 


( 22 ) 


( 23 ) 


respectively. Therefore, in the tangent moduli formulation, the material param- 
eters appear inside the integrals while in the secant moduli formulation the material 
parameters are outside the integrals. 

For an elastic response the formulation in Eqs. (2) or (23) is preferred but 
for an inelastic response the form has not been determined based on first principles. 
The material constants can appear either inside or outside the integrals. The 
correct choice can be determined by comparing analytical predictions for either choice 
with experimental results. This will be the subject of the following sections for 
the three specific constitutive models considered. 


3.3 Discussion of the Three Constitutive Formulations 

The three constitutive theories to be considered are: (1) Walker's, (2) Miller's 

and (3) Krieg, Swearengen, and Rhode's. Walker's Theory, from Ref. 1, is summarized 
in Appendix 3 in differential form and Appendix 4 in integral form. Equation (3.2) 
in Appendix 3 contains rate of change of temperature terms. These terms are a 
result of assuming the material constants n^ and n£ appear outside the integrals 
in Eq. (4.2) of Appendix 4. Miller's equations presented in Appendices 5 and 6 
from Ref. 1 have an assumed temperature dependence based on activation energies. 
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•a 

jff 


Therefore, the material constants are not functions of temperature. An adjustment 
in the theory is then necessary to include rate of change of temperature terms. 

Krieg, Swearengen and Rhode's Theory from Ref. 1 is summarized in Appendices 7 and 
8. Again, this theory contains no rate of change of temperature terms but the 
necessary modification can be made simply by moving the appropriate constants out- 
side the Integral. Tables 1, 2 and 3 present the material parameters for each of 
the three theories used in Ref. 1. 

* 

3. A Rate of Change of Temperature Effects in Walker's Theory 

Walker's Theory will produce results during nonisothermal loading that depend 
on rate of change of temperature terms, in addition to the effects during elastic 
loading. This is reflected by Eqs. (3.2) and (4.2) in Appendices 3 and 4. With 
respect to Walker's theory, as reported in Ref. 1, the differential form has been 
modified to o more accurately reflect the integral form by the addition of the time 
derivative %-j. With respect to the elastic constants. Walker's theory requires 
the coefficient of thermal expansion not to be a function of the strain. In the 
remaining two theories, the coefficient of thermal expansion must be a function 
of the strain if the elastic moduli vary with temperature. 

The constants in Walker's theory (taken from Ref. 1) are presented in Figs. 

4-7 and in Table 1 as functions of temperature. Figures 4-7 are used as an aid in 
determining the proper interpolation functions and in determining derivatives with 
respect to temperature. These material constants exhibit rapid changes in the tem- 
perature range of 1000 F to 1400 F. The most important constant in the rate of 
change of temperature terms is n£ (see Fig. 7 and Eq. (3.2) of Appendix 3). The 
parameter, n 2 , displays a reversal in slope in the range from 1000 to 1400 F. 


The user subroutine HYPELA of the MARC code, from Ref. 1, was modified to in- 
clude rate of change of temperature terms. Figures 8-13 present the results from 
Walker’s theory for the three TMF loops of Figs. 1-3, with and without the rate 
of change of temperature terms. Although the general shape of the loop remains the 
same, the steady state hysteresis loop is approached much more rapidly with the 
rate of change of temperature terms than without. 


3.5 Rate of Change of Temperature Effects in Miller's Theory 

Miller's theory, presented in Appendices 5 and 6 and taken from Ref. 1, was 
modified by making in Eq. (5.2) temperature dependent. The modified theory is 
summarized in Appendix 9. Figures 14-19 illustrate hysteresis loops at different 
tempratures for the values of listed in Table 4. These hysteresis loops should be 
compared to the corresponding experimental loops presented in Figs. 20-25 taken from 
Ref. 1. The remaining constants are taken from Ref. 1, and summarized in Table 2. 

In Miller '8 theory, with the constant now assumed to be a function of tem- 
perature, the subroutine HYPELA in MARC was modified to include rate of change of 
temperature terms. 
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For the three THF loops, Miller's Theory, with and without the rate of 
temperature terms, Figs. 26-31 show almost no change. This is partly due to 
fact that steady state is approached rapidly even without the rate of change 
of temperature terms and partly due to the fact that the material constant does 
not vary appreciably with temperature between 800 and 1200 F. It is difficult 
using Miller's theory to arrive at empirical values of at low temperatures that 
match the experimental results (for example, compare Figs. 14 and 20). Therefore, 
the lower temperature estimates for can vary appreciably from the chosen value. 

3.6 Rate of Change of Temperature Effects in 
Krieg, Swearengen and Rhode’s Theory 

Krieg, Swearengen and Rhode’s theory, Appendices 7 and 8, was modified by moving 
out of the integral in the integral form of the theory. Appendix 10 presents 
a summary of the modified theory, and Table 3 summarizes the constants which were 
taken from Ref. 1. 

These modifications were included in a version of the MARC user subroutine 
HYPELA (see Ref. 1). The modified subroutine was tested using the three thermo- 
mechanical loops. Figures 32-37 present the resulting stress-strain hysteresis 
loops with and without the rate of change of temperature terms. For the two open 
thermomechanical loops, there was little change when the rate of change of tem- 
perature terms were included, compare Figs. 33 and 34 with Figs. 36 and 37, respec- 
tively. For the closed thermomechanical loop, Figs. 32 and 35, the steady state 
loop is approached more rapidly when the rate of change of temperature terms are 
included. 


3 . 7 Summary 

In all succeeding analyses only one viscoplastic theory is used. For these 
analyses Walker's theory was chosen becasue of its ability to represent more phys- 
ical phenomena than Krieg, Swearengen and Rhode's theory. Miller's theory with its 
assumed variation of the material parameters with temperature made it more difficult 
to represent experimental hysteresis loops. A more recent version of Miller’s theory 
Ref. 2 and Appendix 11, removes some of this difficulty but the number of material 
parameters has been substantially increased. 
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4.0 TASK II - TECHNICAL DISCUSSION: TIME INDEPENDENT RESPONSE 


The Task II objective is to incorporate into one candidate theory from Task I, 
a new state variable to accommodate an inelastic instantaneous response. For 
comparative purposes, the three thermomechanical "faithful cycle" hysteresis loops 
analyzed in Task I are used in this task. TVo forms for representing instantaneous 
response will be discussed. The first is based on an instantaneous stress state 
variable. The second uses an inelastic work rate. The last section will present 
results based on the inelastic work rate. 


4.1 Instantaneous Stress State Variable 

A particular formulation for the instantaneous stress state variable has 
been examined for uniaxial stress states. The formulation is summarized in 
Appendix 12. The particular case considered is summarized in Table 5, and includes 
all of the important features of Walker's formulation. The equations under these 
assumptions reduce to 

. .j. _ / lor— XII \ n ___/_ ^ 
cr- 


Cl K J 

(24) 

■ . & 


c “ r 

(25) 

3 «m 2 c-m 3 lcls 

(26) 

ft “n 2 c-n 3 lclft 

(27) 


*1 
. I 


Consider the case of initially loading the material at a specified strain rate 
and take 
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€*c = o, XI ■ H m o' =o,k -y 


From Eqs. (24-27) 



t = 0 


(28) 


(29) 


and the apparent Young's Modulus, E 0 , during Initial loading Is 

JL- ! + -L 

E 0 E m 2 


(30) 


If the material is. new loaded at constant strain rate, y, until steady condi- 
tions result, and then suddenly the strain rate is changed to -y, the material 
will now unload initially by 



(31) 


which can be evaluated in a similar manner from Eqs. (24-27). Then the apparent 
Young's Modulus during fully reversed unloading, E rev » is 


I 


E 


rev 



+ 


1 

2m 2 


(32) 


therefore, the unloading slope is not equal to the initial loading slope in contra- 
diction to experimental results. In order for the two slopes to be nearly equal, 
the constant m 2 must be large. A large value for m 2 should result in the instan- 
taneous stress quickly saturating to its steady value and then having no subse- 
quent effect on the response. This formulation is therefore not acceptable. 


4.2 Inelastic Work Rate Formulation 

A formulation based on the rate of work being done appears to be able to 
represent some of the time independent features appearing in the experimental data. 
The time independent terms can be developed based on a modification of the equations 
of Ref. 3. In Ref. 3, Walker used an expression of the form 
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w - « _ - 


r mox 


• t> & 

" a E 


(33) 


where a Wf>y , b and a are positive constants. As b/a approaches unity, the stress- 
strain curve on initial loading becomes identical to an elastic perfectly plastic 
response. There are limitations with this formulation, two of which are: 


(1) the initial loading slope (i.e.. Young's modulus) is not equal to an 
Initial unloading slope if a is not equal to b, and 


(2) the inelastic strain rate, c, is opposite in sign to the total strain rate, 
e, during Initial unloading (i.e., if a > o then c > o even if e < o) . 


Equation (33) does produce a time independent response similar to classical 
plasticity and therefore modifications of this formulation were developed. In order 
to incorporate a time independent formulation into viscoplastic theories, the 
inelastic strain rate, c^j, is decomposed into two parts as follows 


r..= c®+ c p . 
C 'J C 'J C 'I 


( 34 ) 


c p 

where c^j is the time dependent part, similar to classical creep strains, and c^ 

is the time independent part, similar to classical plastic strains. 

Such a decomposition is being developed by Haisler, Ref. 4, and Bradley, 

Ref. 5. Using Walker's theory, the time dependent inelastic strain is 


fyfi 


(f’if “ij> 


(f s if n ii 

) 

" 1 





( 35 ) 


The time independent inelastic strain rate can be taken as an expression similar to 
Eq. (33) but using work rates rather than strain rates. The total rate at which 
work is being done would replace the total strain rate. The elastic strain rate, 
o/E, would be replaced by a recoverable work rate. The total work rate, W, is 
given by 


" s Vij 


( 36 ) 
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and the nonrecoverable work. 




9 


will be taken as 




ORIGINAL PAGE fS 

OF POOR QUALITY 


( 37 ) 


where < x > ■ \ is the unit ramp function, 

and can be used to represent kinematic hardening, 
strain can be taken to be 


has a zero trace (u^ = o) 
The time independent inelastic 




2<r* 

oo 


< w-k (w-w p )> 


(38) 


where is the time independent yield stress and can be taken to be a function, 
for example, of temperature, nonrecoverable work (to simulate isotropic hardening), 
and if necessary total strain rate. 


Substituting Eq. (38) into Eq. (37) 


W Pj §■ <(l-k) W+kw P > 

2o- 2 

oc 


(39) 


where 3J£ ■ \ S^ - s ij" w ij)* Equation (39) can be solved for 0 P . 

3(l-k)j' 




n 1 1 

. n oo 

W K ! 7 TT- 

- m 


(40) 


where use has been made of Eq. (37). Substituting Eqs. (36) and (40) into Eq. (38) 
gives the time independent inelastic strain 


KSn-fwij) it-kXffjjijj* 

c 'i° 


" k 


(41) 


For classical plasticity without hardening, k is unity and u. . vanishes. Equation 
(4i) becomes indeterminate if 
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J 5 T«r 2 


ORION',?. pr.w k? 

OF PC/'*-; " v / ' 


Therefore, o B can be identified with the classical uniaxial yield stress. Under 
these conditions, Eq. (38) becomes 

3S; j 


which is identical to classical plasticity theory as shown in Appendix 13. 

Appendix 14 presents Walker's theory modified to include these time independent 
strains with no associated hardening. There are two additional material constants, 
k and o,*,, which are functions of temperature. Equation (41) in its most general 
form is capable of representing a wide range of stress-strain behavior. For 
example, taking 

w ii = Hc »i P 
'J 


to represent kinematic hardening, and 


<5.= «o + <r, 


V*/q 

1+qV 


I + q. fin 


where H is the kinematic hardening slope, 
o 0 is the Initial "yield" stress, 

and 02 can be used to represent isotropic hardening, 
q^ can be used to represent strain rate effects, and 
e 0 is the strain rate appropriate to o Q . 

A uniaxial loading simulation was performed using the equations in Appendix 14, 
but with Eqs. (44 and 45) together with the data in Table 6. Figure 38 presents the 
resulting stress strain hysteresis loop which is similar to the experimental loops 
at 800 deg F presented in Fig. 22, and Fig. 39 presents the inelastic strain 
as a function of total strain hysteresis loop. 
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4 . 3 Time Independent Hysteresis Loops 

The time independent response equations can be simplified by using the rest 
stress, ft (equilibrium stress), to represent kinematic work hardening in the time 
dependent response and proportionally varying the equivalent yield stress, o*, in 
the same manner as the drag stress, K. Appendix 13 summarizes the governing equa- 
tions. This then leaves two additional constants to be evaluated for the theory 
o,*, (at R equal to zero) and k. The constant k has been taken to be zero for all 
temperatures while o M was set to infinity for 1200 F to 1800 F. At 1000 F and 800 F, 
the constant o w was determined and the constants n 2 and n 3 were changed to more 
closely match the experimental results. In addition, the constants and 1/n at 
1400 F were changed to more closely match the experimental results. The 1400 F 
change only effects strain rate results higher than those presented. Table 7 
summarizes the new values for each of the constants. Figures 40-47 show the 
resulting hysteresis loops. At 1000 F, the loops are significantly more accurate 
with the new values of n 2 and n^» compare Figs. 20-25 with Figs. 40-47, than the 
predictions of Ref. 1. The effect on the predicted stress when the time independent 
terms are included will be discussed in the next section. 

The time dependent and time independent inelastic strains can be separated 
experimentally by performing creep and relaxation tests beginning various points 
on the steady state stress strain hysteresis loops. 
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5.0 TASK III - TECHNICAL DISCUSSION: COMBINED TASK I AND TASK II EFFECTS 


The Task III objective is to evaluate the theories developed to determine the 
combined improvement in predicting stresses and strains. The analytical results of 
this task are compared with those of Tasks I and II and Contract NAS3-22C55 uting 
the same thermomechanical hysteresis loops and geometry models used in Task I. 

The subroutine HYPELA was modified to include the time independent terms as a 
part of Task II. The material constants in their final form appear in Table 7, 
the initial form of the constants appear in Table 1. Figures 48-53 show the effect 
of changing constants and including time independent terms and rate of change of 
tcuperature terms for the closed symmetric TMF cycle. The other two TMF cycles 
were only minimally effected (see Figs. 54 and 55). Figures 48-53 show a definite 
improvement in accuracy primarily due to the change in the low temperature material 
properties . Adding the time independent terms and the nonisothermal terms has only 
a minimal effect. 

A slight improvement in the predicted response of the open nonsymmetric TMF 
cycle was attained by setting the material constant o« at 1200 F to 60 ksi with no 
appreciable change in the other two TMF loops, as illustrated in Figs. 56 to 58. 
Figure 59, taken from Ref. 1, shows the predicted stress strain response for the 
open nonsymmetric TMF cycle using the material parameters from Table 1. Comparing 
Figs. 58 and 59, there is a considerable improvement in the predicted response when 
compared to the experimental response (Fig. 3). This is primarily due to the change 
in low temperature material constants. 

Summarizing the combined effects of Tasks I and II showed that accurately 
modeling the low temperature hysteresis loops had the most effect, followed by the 
instantaneous inelastic response. The rate of change of temperature terms studied 
under Task I produced changes in the thermal ratcheting. Thermal ratcheting seemed 
to occur more often when there were sudden changes in the slope of the stress strain 
curves . 


18 



R8 3-956077-1 


6.0 TASK IV - TECHNICAL DISCUSSION: IMPROVEMENTS IN 

COMPUTATIONAL EFFICIENCY 


The Task IV objective is to improve the computational efficiency by: 

(1) incorporating an array containing the state variables into subroutine 
HYPELA of the MARC program to speed up computation time, 

(2) further reducing computation time by modifying the MARC program so that the 
stiffness matrix is assembled once only during a calculation, with no 
further reassembly during succeeding increments. Further, more 
efficient programming of the computer models developed in Tasks I and II 
and Contract NAS3-22055 (Ref. 1) is completed, and 

(3) performing test cases of sufficient complexity to adequately demonstrate 
the improvements in computational efficiency. 

The first part of this section presents various numerical techniques for 
integrating the differential equations for the modified Walker theory summarized 
in Appendix 14. The second part summarizes the computational savings that results 
when various modifications were made to the MARC finite element code and the user 
subroutine HYPELA of the MARC code. Reference 6 presents a complete discussion for 
the modifications to the MARC code necessary to incorporate the equations of 
Appendix 14 in an efficient manner. 


6.1 Numerical Integration of the Modified Walker's Theory 


Three methods for integrating the constitutive equations were examined: 

(1) a forward difference with error estimates for revising the time step, (2) a 
backward difference, and (3) integrating an integral form of the equations. 


The forward difference Integration, similar to Ref. 
estimate, E, given by 


E s AR + 



7, is based on an error 


( 46 ) 


where 





£- 


,t&' >t 
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and 


A V i as ij A3 ij 


A of a quantity is the change in the quantity over time step At. 
If the error estimate is too large 


«>«. 


( 47 ) 


then the time step is halved and the step is repeated. If the error estimate is 
too small 


« < € 


2 


( 48 ) 


then the tim°. step is doubled for the next integration step. 

Figure 60 shows the resulting stress-strain loop for the symmetric closed TMF 
cycle using a 200 step per cycle forward difference without the error estimate. 
This should be compared with Fig. 61 where 


and 


«,s I0* 4 

V'°' s 


For this case, 195 steps per; cycle resulted with an increase in computation time of 
0.4 percent. In Fig. 62, the allowable "errors" have been increased ten times and 
although the response demonstrates some instability, the stress-strain loop is still 
approximately correct. In Fig. 62, only 38 steps per cycle were required and a 
decrease in computation time of 75 percent resulted. This is obviously a viable 
method for determining the time step magnitude automatically. 

A backward difference algorithm was developed which is summarized in Appendix 
15. Figures 63-66 display the resulting stress-strain loop for the closed 
symmetric cycle for various numbers of increments per cycle. At 16 Increments per 
cycle, Fig. 65, the loop is marginally accurate when compared to Fig. 60. This 
represents a savings of approximately 75 percent 6ince about 2.5 iterations to solve 
the resulting nonlinear simultaneous equations per loading increment were required. 
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Although the forward difference worked adequately for all tesc cases considered, 
the backward difference convergence was slow in cases where the term n* was not 
equal to zero and the strain rate was small (i.e., on the order of 10“*/sec). A 
quadratic Newton’s method was also used to solve for LG and but there was no 

benefit over a linear Newton's method. Therefore, in these cases, a sufficiently 
small step size was required. 

The equations for the modified Walker's theory were recast in integral form, 
see Appendix 16. This form was also tested and the results were similar to the 
results for the backward difference. 

Summarizing, three methods for integrating the equations were examined They 
included; (1) forward differences using error estimates to revise the time step, 

(2) backward differences, and (3) recasting the equations ?.n integral form and 
then numerically integrating. Tests of these methods were .’formed on the one 
dimensional form of the equations. The forward differences with error estimates 
saved up to 75 percent of the computational time. The backward difference integra- 
tion was stable for large time steps and converged rapidly for typical thermal cycle 
problems, resulting in a savings of computational tim-s of about 75 percent. Un- 
fortunately, at low strain rates and high temperatures, where thermal recovery can 
be important, convergence was sometimes unacceptably slow. Results for the integral 
form indicated that the savings and problems are similar to the backward differences. 
The convergence problems make the backward difference and integral form integration 
methods unreliable and therefore only the forward difference integration with error 
estimates was used in the MARC code modifications. 


6.2 Computational Efficiency Improvements 

Three approaches were investigated to improve the computational efficiency. 

One of these is the numerical integration discussed above. The second included the 
approach of storing the state variables in an array in HYPELA. This proves to be 
advantageous only if all of the 16 state variables for each integration point of 
each element can ba stored in core. Therefore, there is a savings only on relatively 
small models. In the third approach, the MARC code was modified to save the inverted 
stiffness matrix (i.e., use back-substitution), thereby saving the CPU time necessary 
to solve the simultaneous equations. 

In the subrouinte HYPELA of the MARC code, the constitutive equations were inte- 
grated using a forward difference for the time derivative. Each increment in the 
MARC code is divided into a number of subincrements. The number of subincrements 
is stored as the variable NSPLIT. For an ordinary forward difference NSPLIT is a 
constant. For a forward difference integration with error checking, a self-adaptive 
scheme, the variable NSFLIT is charged according to whether the estimated error is 
larger or smaller than the allowable b^ jnds (Eqs. (46-48)). 


"TT 
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The MARC code could then be run in four modes depending on whether NSPLIT 
is constant (ordinary forward difference) or variable (self-adaptive forward 
difference) and whether matrix inversion is performed or back substitution is used. 
IVo different test cases were used in the tests. The first was a two element by 
two element uniaxial axisymmetric mesh, with the loading given by the open non- 
symmetric loop in Fig. 3. The second test case was identical except a seven by 
seven element mesh was used. The first case was run for 48 loading increments and 
the second for 10 increments. Table 8 summarizes the resulting UNIVAC 1110 CPU 
time for five different combinations. Most of the savings resulted when a variable 
time step was used. Using matrix inversion produced a negligible increase in the 
CPU time when compared to back substitution. In the test runs using a constant 
time step the value of NSPLIT was taken to be the maximum value attained in the 
equivalent self-adaptive run. Figure 67 shows NSPLIT at the end of each MARC 
Increment using the 2x2 finite element mesh and self adaptive integration, for the 
first complete cycle of the open nonsymnetric loop. A total of 219 subincrements 
were required. If NSPLIT was set at 16 subincrements per increment then a variable 
time step would produce a 72 percent savings. 

In Fig. 67, at increments 10 to 15 and increments 38 to 45, NSPLIT is large. 

At these increments the larger changes in mechanical and thermal loading occurred. 
In increments 20 to 36 NSPLIT is equal to one which occurs during the hold time for 
the temperature and strain (i.e., the time during which temperature and strain are 
constant) . 

Summarizing, most of the CPU savings was attained by using variable step 
integration while utilizing back substitution produced a relatively small savings. 
Indications are that large problems would have to be run (1000 degrees of freedom) 
before the equation solution time becomes large. 
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7.0 CONCLUSIONS 


During the course of this investigation several significant conclusions were 
drawn. They are: 

1. The variation of material parameters with temperature must be carefully 
considered. This variation is best illustrated by examining the variation of the 
elastic moduli with temperature. If instantaneous, or tangent, moduli are used, 

a perfect elastic reaction where the material undergoes changes in temperature can- 
not be simulated. If secant moduli are used, then elastic nonisothermal reactions 
can be simulated. Similar considerations must be examined for other material param- 
eters in any constitutive model. The use of secant material parameters in time 
dependent reactions introduces terms containing rate of change of temperature. 

2. Including rate of change of temperature terms in the equilibrium stress 
formulation can have a significant effect on the stress or strain, ratcheting 
during cyclic thermomechanical fatigue loading. The thermal ratchetting was also 
found to be sensitive to the material parameters used. 

3. An instantaneous (i.e., plastic) response term has been introduced which 
can be shown to reduce to classical perfect plasticity. This term allows more 
accurate reproduction of low temperature material response. 

4. The importance of accurately evaluating the material parameters is 
illustrated by comparing the poor comparisons of the end of Task I, Fig. 9, with 
the improved results at the end of Task III, Fig. 50. 

5. A sixty percent reduction, Table 8, in CPU time can be attained by 
integrating the constitutive equations using a variable time step when compared 
to a fixed time step. 
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TABLE 1 



Material 

Constants 

for Walker* 

s Theory 

(Ref. 1) 


Material 

Constant 

982°C 
1800 °F 

871°C 
1600 °F 

760 °C 
1400°F 

648°C 
1200 °F 

537°C 
1000 °F 

427°C 

800°F 

A 

11.5E6 

15.4E6 

17.8E6 

18.1E6 

17.2E6 

17.8E6 

V 

4.9E6 

6.9E6 

8.4E6 

9.0E6 

9.0E6 

9.8E6 

K i 

59292 

91505 

251886 

95631 

75631 

50931 

K 7 

0 

0 

0 

0 

0 

0 


\ 

-l 

n 

.233 

.195 

.244 

.079 

.059 

.059 

i 

m 

1.16 

1.16 

1.16 

1.16 

1.16 

1.16 

1 

i 

n l 

0 

0 

0 

0 

0 

0 

< -V 

-L ( t 

n 2 

1.0E6 

5.0E6 

2.0E7 

1.5E7 

6.0E7 

30.0E7 

i 

i 

] 

n 3 

312 

673 

1179 

781 

1000 

8000 

vn 1 

► •• 

n 4 

0 

0 

0 

0 

0 

0 

: 1 
; .r 

n 5 

0 

0 

0 

0 

0 

0 

, ! 

/' i 

n 6 

2.73E-3 

8.98E-4 

0 

0 

0 

0 

• i 

i 

n 7 

0 

0 

0 

0 

0 

0 

* ! 

8 

-1200 

-1434 

-2000 

-2000 

0 

0 


i 



f 
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TABLE 2 

Material Constants for Miller's Theory (Ref. 1) 


B 


H, 


H„ 


2 

Q* 

T 

m 

k 

e' 

e' 


8000 

1.598 

1.0293E14 

1.0E7 

9.305E-4 

100 

50000 

5.9425E-12 
104600 
1588 °K 


1.9859 
exp 
exp 


for T * . 6T, 


Q*/kT | 


m 


for T < .6T 


m 


where T is the temperature in degrees Kelvin. The Lame constants X and y 
are as given in Table 1. 
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TABLE 3 

Material Constants for Krieg, Swearengen and Rhode's Theory (Ref. 


Material 

982 °C 

871°C 

760°C 

648°C 

537°C 

Constant 

1800°F 

1600 °F 

1400°F 

1200°F 

1000 °F 

A 

11.5E6 

15.4E6 

17.8E6 

18.1E6 

17.2E6 

V 

4.9E6 

6.9E6 

8.4E6 

9.0E6 

9.0E6 

K 

0 

59292 

91505 

251886 

95631 

75631 

n 

.223 

.195 

.244 

.079 

.059 

A i 

1.0E6 

5.0E6 

2.0E7 

1.5E7 

6.0E7 

A 2 

243 

14.96 

1.54 

.66 

1.79E-3 

A 3 

1.0E-12 

1-0E-12 

1.0E-12 

1.0E-12 

1.0E-12 

\ 

0 

0 

0 

0 

0 

A 5 

0 

0 

0 

0 

0 


427°C 
800 °F 

17.8E6 

9.8E6 

50931 

.059 

30.0E7 

.59 

1.0E-12 

0 


0 
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b 


TABLE 4 


Constant 

in Miller's Theory 

Temperature 

H 1 

deg F__ 

£Si 

800 

10 x 10 6 

1000 

10 x 10 6 

1200 

8 x 10 6 

1400 

10 x 10 6 

1600 

6 x 10^ 

1800 

3 x 10 6 


r 



£■ 
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TABLE 5 

Parameters for Instantaneous Stress State Examination 


Parameter 


K, 


n l* “l 


V m 2 


n 3 , m 3 


n 4-6’ m 4-6 


Value 

nonzero constant 
nonzero constant 
nonzero constant 
0 
0 

nonzero constant 
nonzero constants 


0 
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TABLE 6 

Constants in Time Independent Simulations 
k * 0.3 


f 

! 



K = 00 

e 0 . 006 

mar. ■ 

q x - .004 

q 2 « 0 

e = .001 

o 

o * 70 ksi 

o 

a 1 - 30 ksi 

H * 7.5 x 10^ psi 






Material 

982°C 

871°C 

760°C 

648°C 

537°C 

427°C 


Constant 

1800 °F 

1600°F 

1400°F 

1200°F 

1000°F 

800 e F 


X 

11.5E6 

15.4E6 

17.8E6 

18.1E6 

17.2E6 

17.8E6 


y 

4.9E6 

6.9E6 

8.4E6 

9.0E6 

9.0E6 

9.8E6 

) 

1 

K l 

59292 

91505 

110696 

95631 

75631 

50931 

I 

K 2 

0 

0 

0 

0 

0 

0 

t 

, 1 

> J 

-1 

n 

.233 

.195 

.1497 

.079 

.059 

.059 


in 

1.16 

1.16 

1.16 

1.16 

1.16 

1.16 

V 

n l 

0 

0 

0 

0 

0 

0 

v , 

-;'V 

n 2 

1.0E6 

5.0E6 

2.0E7 

1.5E7 

1.9E7 

1.0E7 


n 3 

312 

673 

1179 

781 

320 

250 


n 4 

0 

0 

0 

0 

0 

0 


"S 

0 

0 

0 

0 

0 

0 

- 

1 

n ft 

I 6 

2.73E-3 

8.98E-4 

0 

0 

0 

0 

» 1 

n 7 

0 

0 

0 

0 

0 

0 

% 

\ o 

1 ® 

-1200 

-1434 

-2000 

-2000 

0 

0 

*h*> 

V 

*-* 

o. at R • 0 

00 

00 

00 

00 

48000 

48000 

* 

Mt-- 

k 

0 

0 

0 

0 

0 

0 

» 1 

i 

* i 

C V 







V 

V*j' 








\ 
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TABLE 8 

UNIVAC CPU Times 
TMF1 - MARC Increments 0-10 


7X7 Elements 


Matrix Inversion 

NSPLIT 

CPU Time (Sec) 

Percent Savings 

Yes 

Constant 

335 

0 

Yes 

Variable 

141 

58 

No 

Variable 

131 

61 


TMF1 - MARC 

Increments 0-48 



2X2 

Elements 


Matrix Inversion 

NSPLIT 

CPU Time (Sec) 

Percent Savings 

Yes 

Constant 

233 

0 

No 

Variable 

90 

61 




Figure 1. TMF Cycle No. 1: Closed Symmetric Cycle 



700 



Figure 2. TMF Cycle No. 2: Open Symmetric Cycle 








Figure 3. TMF Cycle No. 3: Open Nonsymmetric Cycle 
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Figure 4. Evaluated Power Law Exponent, n, has a Step Change Occuring 
Between 1400 F and 1200 F 
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j Figure 5. Evaluated Thermal Recovery Equilibrium Stress Parameter, n 3 , has a 

| Reversal in Slope Occurring Between 1000 F and 1400 F 
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Figure 7. Constant Back Stress Parameter, K<, has a Major Discontinuity 
Occurring Between 1200F and 1600F 


82 - 4 - 57-12 




ORIGINAL PAGE 
OF POOR QUALITY 


A\ 




< o 

a: g 

co £ 


ed$N 6 9=!S>< l 
08* SS3diS 


82 - 4 - 57—22 


Rate of Change of Temperature Terms 










Rale of Change of Temperature Terms 



Rate of Change of Temperature Terms 





ORIGINAL PAGE IS 
OF POOR QUALITY 



Rate of Change of Temperature Terms 








180 


ORIGINAL PAGE 19 
OF POOR QUALITY 



• 3 — 2 — 87—2 


Figure 15. Hysteresis Loop for Hastelloy-X Using Miller’s Theory at 538 C for a Strain Rate 
of 6.67 x 1<T 3 





Figure 17. Hysteresis Lc.-p for Hastelloy X Using Miller’s Theory at 760 C for a Strain Rate 
of 3.70 x 10* 4 / sec 
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Figure 22. Experimental Steady State Hysteresis Loops for Hastelloy-X at 649°C 
(1200°F) 
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Figure 32. Predicted Response using Krteg, Swearengen and Rohde's Theory for Closed 
Symmetric TMF Cycle Without Rate of Change of Temperature Terms 





Figure 33. Predicted Response Using Krieg, Swearengen and Rohde’s Theory for Open 
Symmetric TMF Cycle Without Rate of Change of Temperature Terms 
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Figure 35. Predicted Response Using Krieg, Swearengen and Rohde’s Theory for Closed 
Symmetric TMF Cycle With Rate of Change of Temperature Terms 
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Figure 36. Predicted Response Using Krieg, Swearengen and Rohde’s Theory for Open 
Symmetric TMF Cycle With Rate of Change of Temperature Terms 





Figure 37. Predicted Response Using Krieg, Swearengen and Rohde’s Theory for Open 
Nonsymmetric TMF Cycle With Rate of Change of Temperature Terms 
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Figure 39. Typical Inelastic Strain Response Using Time Independent Formulation 
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Figure 40. Modified Walker’s Theory at 538 C and 0.0000667/sec Strain Rate 






Figure 41. Modified Walker’s Theory a! 536 C and 0.02/sec Strain Rate 










Figure 42. Original Walker’s Theory at 538 C and 0.000667/i »c Strain Rate 
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Figure 43. Original Walker’s Theory at 538 C and 0.02/sec Strain Rate 




ORIGINAL PAGE IS 
OF POOR QUALITY 

r d ? * g ! o 


Z c 

2 I 


z « z < 


BdW 6'5 = !S>( l 
m 'SS3H1S 


85 - 4 — 58—6 


Figure 44. Original Walker's Theory at 538 C and 0.000066«>sec Strain Rate 





Figure 45. Original Walker’s Theory at 538 C and 0.02/sec Strain Rate 
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Figure 46. Modified Walker’s 1 neory at 427 C and 0.02/sec Strain Rate 














Figure 49. Hysteresis Loop for Closed Symmetric Cycle 
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Figure 50. Hysteresis Loop for Closed Symmetric Cycle 
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Figure 51. Hysteresis Loop for Closed Symmetric Cycle 
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Figure 52. Hysteresis Loop for Close > Symmetric Cycle 
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Figure 54. Hysteresis Loop for Open Symmetric Cycle 
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Figure 57. Predicted Stress-Strain Response of TMF Cycle No. 2 for (649 C) Equal to 60 ksi 
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Figure 58. Predicted Stress-Strain Response of IMF Cycle No. 3 for (849 C) Equs 1 to 60 ksi 






Figure 60 . Forward Difference Integration 
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Figure 61. Forward Difference Integration 
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Figure 63. Backward Difference Integration 
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Figure 64. Backward Difference Integration 












Figure 66. Backward Difference Integration 
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Figure 67. NSPLIT at End of increment for Open Non-Symmetric TMF Cycle 
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Constraints on Coefficient of Thermal Expansion 
Using Tangent Moduli 


For an elastic response 


^ii = ?M«ij + ^<KK 8 »r al3X+2AA> ^ 8i ' 


( 1 . 1 ) 


But is not history dependent, therefore. 


°ij = Fjj (*»T) 


( 1 . 2 ) 


• dF ii • *il * 

a i) = d€ Kt € *t + dT T 


(1.3) 


Comparing Eqs. (1.1) and (1.3) 


bi * 2 ^ 8 «* ^if + X8 'i 8 ** 


(1.4) 


= - a(3X+2^.) 


(1.5) 


The mixed partial derivatives of Eqs. (1.4) and (1.5) must be equal or 

4\ . 2 -£- 8 ik 8, ( ♦ iis, 8 k( . - A^hllA 8 

de Kf dT dT ,H I* dT ,J k * de H1 11 

Contracting Eq. (1.6) on the indices i and j 


( 1 . 6 ) 


l dk a d(aK) 

Ta 7 8 « + 7 ^ 7 =° 


(1.7) 
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where 
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K = X+ ~-fi 


3(1-2 v) 


( 1 . 8 ) 


If the material is linear elastic, the bulk modulus is not a function of strain 
and Eq. (1.7) becomes 


+ JSL s0 

3K dT d« kl 


(1.9) 


Multiplying Eq. (1.6) by the permeation tensor e results in 

Pij 


—■ Epmn= o 


(1.10) 


for all p, m, and n and therefore 


(1.11) 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


Constraints on Coefficient of Thermal Expansion 
for Mixed Formulation 


For an elastic response 


t V (X+ T^ ) {‘kk~ 3 /' t(0 !? d *} V^ir-fV** 


( 2 . 1 ) 


The integral in Eq. (2.1) can be written as 

r* -.1 f\ <Jt i 


r »tn 01 *(- * ^ o~Kk 

J 0 ^ ae 3 tf kk" 3(3X+2^) 


(2.2) 


Taking the derivative with respect of time of Eq. (2.1) and substituting Eq. (2.2) 
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(2.3) 


But a is not history dependent or 


<T 'j = F 'i 


(2.4) 


Taking the derivative with respect to time 
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ai i‘a« hi € ki + aT T 


(2.5) 


Comparing Eqs. (2.3) and (2.5) 
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j (2.6) 
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Appendix 5. Miller’s Theory (Differential Form) 





is the initial value of K, and T is the lemperalure in *i 


( 




b & 


ORIGINAL PAGc £3 
OF POOR QUALITY 


— ScP 
£ * 


2 cT 

** II 
« rf 

£> 5. 


jy o 
m x: 
5C I 
CM f* 

< 5^ 
^ sc 





\n 


cp 

r- 


C\) 

r*- 


CO 



ORIGINAL PAGE 13 
OF POOR QUALITY 


o 

o 

3 

15 

> 


2 


0 
* 
0> 

1 

5 




I 


*1 



i 



depend on temperature 





ORIGINAL PAGE 13 
OF POOR QUALITY 



Material constants, n * Hj , A, , A? t Cg-.O*, k are independent of temperature. 

Material constants A, M ,1 <o .0 .H, depend on temperature, where K Q is the initial value of K. and T is the temperature i 


Appendix 10. Modified Krieg, Swearengen and Rohde’s Tneory (Differential Form) 
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Material constants: K^i,n ,A,,A 2( A3 ,A 4 ,A 5 ,K Q depend on temperature, where K 0 is the initial value of K. 
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Appendix 11. Updated Miller’s Theory - Dllferenllel Form 
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Appendix 11. Updated Miller’s Theory — Differential Form (Cont) 
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Appendix 12. Walker’s Theory — One Dimensional Differential Form including Instantaneous 

Stress State Variable 
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Parameter 8 is: 0 if instantaneous stress is inelastic 

1 if instantaneous stress is elastic (previous formulation) 
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Plastic Strain Rates in Classical Plasticity 



Consider an elastic perfectly plastic material, 
is given by the normal to the yield surface 


the plastic strain rate, 




whare 




1 s « 

■ 2 °kk S lj 


is the second invariant of the deviatoric stress tensor 
, and a is the uniaxial yield stress. Then 





where X is a scalar that must be determined. 


The plastic work rate is 


w p=: 0'j J *Cjj P = X s ij s ij = 2 J 2 x 




and therefore 


W p _ 3w P 
'2J 2 ‘ 2 a y 2 


Substituting Eq. (13.4) into Eq. (13.2) 



(13.1) 


(13.2) 


(13.3) 


(13.4) 


(13.5) 


which is identical to the result presented in the text. 
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Appendix 15. Backward Difference Integration Algorithm 
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Appendix 15. Backward Diffaranca Integration Algorithm (Cont) 
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Appendix 16. Modified Walker’s Theory — Integral Form 
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